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ABSTRACT 

The pioneering study by Skumanich (1972) showed that the rotational velocity of 
G-type Main-Sequence (MS) stars decreases with stellar age according to (vsini) ex 
^-1/2^ This relationship is consistent with simple theories of angular momentum loss 
from rotating stars, where an ionized wind is coupled to the star by a magnetic field. 
The present study introduces a new approach to the study of stellar rotational braking 
in unsaturated F and G type stars limited in age and mass, connecting angular mo- 
mentum loss by magnetic stellar wind with Tsallis nonextensive statistical mechanics. 
As a result, we show that the rotation-age relationship can be well reproduced using 
a nonextensive approach from Tsallis nonextensive models. Here, the index q, which 
is related to the degree of nonextensivity, can be associated to the dynamo process 
and to magnetic field geometry, offering relevant information on the level of stellar 
magnetic braking for F- and G-type Main-Sequence stars. 
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1 INTRODUCTION 

The relationship between stellar rotation and age is an ex- 
citing topic in astrophysics posing, a number of unanswered 
questions, including: 1) How is the stellar rotation-age rela- 
tionship dependent on mass? 2) How does rotation influence 
the stellar dynamo, as well as chromospheric and coronal 
heating mechanisms? 3) How do initial star formation sys- 
tems afi'ect rotational evolution? 4) How do stellar angular 
momentum and stellar activity evolve over time? 

|Skumanich| ( |1972[ )'s pioneering study showed that v 
sini of G-type Main-Sequence (MS) stars for Hyades and 
Pleiades measured by Kraft ([1967) were consistent with 
{vsini) oc t~^", where t is the stellar age, v is the equato- 
rial rotational velocity of the star and i the inclination angle 
of the rotational axis to the line of sight. This relationship 
is in line with simple theories of angular momentum loss 
from rotating stars, where an ionized wind is coupled to the 
star by a magnetic field ( Schatzman|1962 l. Indeed, angular 
momentum loss due to stellar winds is generally believed 
to be responsible for the Skumanich relationship, but the 
exact dependence of rotation on age has not been fully es- 
tablish ed ([Kawaler 1988, K rishnamurthi et aZ.|1997[ ). |Barry| 



et al. (19871 and Sodcrblom et al. (19911 reported similar 



qualitative results for solar-type stars, but with power-law 
exponents ranging from —1/2 (corresponding to the Sku- 
manich relation) to —4/3. [Pace fc Pasquini^ (2004 1 ) claimed 
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that these power-laws do not fit the age-activity-rotation 
of G dwarf stars in open clusters. According to these au- 
thors, a t~^" law is more consistent with the observations. 
Despite the differences in rotation-activity-age relationships 
obtained to date, most of the data strongly suggest that 
this relationship is deterministic and not merely a statisti- 
cal artifact. Indeed, rotation data for stars at different evo- 
lutionary stages and environments, in particular for stars in 
open clusters, have confirmed the deterministic nature (e.g.: 
(jBarnes fc Kim||2010| [Rcincrs and Mohanty||2012[|Barnes| 
(2003), M amajek fc Hillenbrand| ( |2008[ ), [James et gZ | ( |2010[ ) 
and |Mcibom et fflZ.[ ( |2011a|b| ). 



In this study we present a new approach to the study of 
stellar rotational braking in F and G type stars, connecting 



angular momentum loss by magnetic stellar wind (Kawaler 



1988 Chaboyer et aZ.[[1995[) with Tsallis' nonextensive sta- 



tistical mechanics (Tsallis [19881. In Section 2, we revisit 



parametric models for angular momentum loss by magnetic 
stellar wind, with an emphasis on a modified Kawaler model. 
Section 3 describes our nonextensive approach, in Section 4 
we describe the observational data, and Section 5 compares 
our model with v sini measurements for a large sample of F- 
and G-type single main-sequence stars. Finally, conclusions 
are presented in Section 5. 
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2 REVISITING ANGULAR MOMENTUM 
LOSS BY MAGNETIC STELLAR WIND 

Following the theory of angular momentum loss elaborated 
by|Mestel||T968lp84l ) and |Mestel fc Spruit HT987t , |Kawaler| 
1 1988 1 developed a general parametrization showing that the 
rate of angular momentum loss is proportional to fj^^'*"'^"; 
where il is angular velocity, a denotes the dynamo relation- 
ship and A^ is a measurement of magnetic field geometry. 
The problem with angular momentum loss due to magnetic 
stellar wind is its dependence on these two parameters, a 



and A'^. Kawaler ( 1988 1 simplified the dynamo relationship 



by assuming that the total magnetic field is proportional to 
the rotation rate to some power, denoted by Bo oc f2", where 
a varies between 1 and 2 for the unsaturated domain and 
for the saturated domain. As demonstrated by |Chaboyer| 
\et al.\ ( |1995[ ), the A^ parameter exhibits a wide range, vary- 
ing from to 2, where represents dipole field geometry 
and 2 a radial field. Moreover, these authors have shown 
that angular momentum loss from high rotation stars can- 
not be described by Kawaler's relationship. The same au- 
thors and subsequent studies ( [Reiners and Mohanty 2012 1 
have also demonstrated that magnetic fields at high rota- 
tion rates must saturate, following the form dJ/df oc Q 



In addition, studies by Matt and Pudritz (Matt & Pudritz 
|2005[ |2008a|b | have explored the role of stellar winds to the 
spin-up and spin-down torques expected to arise from the 
magnetic interaction between a slowly rotating pre-main- 
sequence star and its accretion disk, demonstrating that if 
stellar mass outflow rates are substantially greater than the 
solar rate, the stellar winds can carry off orders of magni- 
tude with more angular momentum than can be transferred 
to the disk. 

In terms of parametrization, there are two parameter- 
ized forms known as modified Kawaler's relationships: 
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for fl < UJsat- 

Here, u^at is the threshold angular velocity beyond 
which saturation occurs, K^ is a calibration constant and 
R, M and M represent the stellar radius, mass and mass- 
loss rate in units of 10~^*Moyr~'^, respectively. In this re- 
gard, LOsat completes the triplet-order parameter. In fact, 
this saturation threshold is an important ingredient for the- 
oretical models of stellar angular momentum evolution. For 
instance, Krishnamurthi et al. { 19971 and Sills, Pinsonneault 



and Terndrup (20001 suggest that the saturation thresh- 



old depends on stellar mass and is inversely proportional 
to the global convective overturn timescale, which is related 
to the Rossby number, an indicator of the strength of dy- 



namo magnetic activity. By contrast, [Reiners and Mohanty | 
(2012J argue that Usat does not depend on mass. According 



to Pizzolato et al. ( 2003 1 , the saturation phenomenon is also 
based on magnetic proxies and its level can be inferred, irre- 
spective of angular momentum evolution arguments, which, 
do not require a Rossby number for either regime. These 
authors also argue that the X-ray emission level from satu- 
rated stars depends only on bolometric luminosity, whereas 
for non-saturated stars the emission level depends only on 
the rotational period. This implies a direct mass-dependent 
UJsat (jKr ishnamurthi et a/.||l997[ ). Since this timescale does 
not change much in our model, we neglect its effect for a 



given mass (Barnes and Sofia 1996 Krishnamurthi et al. 
1997[ [Pizzolato et fflZ.|2003| ) 



In the present study we assume that since the moment 
of inertia / and stellar radius R changes slowly during MS 
evolution, the angular momentum loss law can be simplified 
by considering the limit of AJ /At for constants / and R, that 
is, this loss law is fully specified by the rotational velocity. 



while the star is braked by the stellar wind (Bouvier, Fores- 
tini fc Allain|1997l. We also consider that in the absence of 



angular momentum loss, equatorial rotational velocity v can 
be determined by the simple condition of angular momen- 
tum conservation, denoted by 



J =4: 



(3) 



It is important to underline that in our model we con- 
sider a solid-body rotation, although this assumption can 
directly affect spin down effects for F-type stars. For a more 
complete model, it is necessary to take in consideration other 



physical aspects such as core-envelope decoupling ( Cameron 
|fc Li Jianke||1993| |Allain|[T998| , as well as the evolution 



ary behavior of the moment of inertia( Wolff fc Simon|1997 
Zorec fc Royer|2012[ ) 



Two pairs of equations were combined to obtain the 
time dependence of v. First, from eqs. M and ([3| we have 

{t-toY 



){t) — vo exp 



with 



T\ 



(to ^ i < tsat) 



(4) 



r,^[mu.t:r'Y" (5) 

for the saturated domain, where /(A) is denoted by 



/(A) = 



R, 



M, 



10-14 



l-2iV/3 



(6) 



where / denotes the moment of inertia of the star and A 
denotes a parameter set given by A = •^ _R, M, M, A'^ >. 

On the other hand, by combining eqs. ([2| and ([3| we 
obtain the equation 



j{t) = Vsat 



it - tsat) 



-3/4aJV 



{t ^ tsat) 



with 



4aA'' , , iaN/a 

—^.fW^sat 



(7) 



(8) 



In this respect, the ri and r are the MS spin-down 
timescales in the saturated and unsaturated regimes, respec- 
tively. The other parameters, to, vq and tsat, are the age at 
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Table 1. Best parameter values of our unsaturated model using 
eq. ||14|l. The values of A'^ using eq. l|15[l also are shown. 



Stars 


(M/Mq) 


Q 


Na^2- 


-Na=l 


F0-F5 


1.36 


l.SOitO.l 


0.30- 


-0.60 


F6-F9 


1.22 


1.96±0.2 


0.36- 


-0.72 


G0-G5 


1.11 


2.18±0.2 


0.44- 


-0.89 


G6-G9 


0.98 


4.38±0.6 


1.27- 


-2.53 


AUF 




1.91±0.1 


0.34- 


-0.68 


AUG 




2.18±0.4 


0.44- 


-0.89 



which the star arrives on the MS, the rotational velocity at 
that time and subsequent age at which the star slows down 



to values below the critical velocity Vsat, respectively (Rein- 



ers and Mohanty 20121. Time teat can be determined by 
setting v(t) = Vsat in eq. Q. Thus, we have 



tsat =to+Tl In 



Vsat 



(9) 



Eq. Q reveals that in the saturated regime the star 
slows down exponentially from time to up to velocity Vsat 
in time tsat- Eq. ([7| shows that from this time onward, the 
spin-down rate decreases to a power law given by —S/AaN. 
Thus, the star remains within a factor of a few Vsat for the 
remainder of its MS lifetime. 



3 NONEXTENSIVE APPROACH TO THE 
STELLAR ROTATION-AGE RELATION 

This section details a nonextensive approach for the modi- 
fied Kawaler model divided into two classes: saturated and 
/3-saturated models. Both have the same unsaturated model, 
but with different saturated timescale tsat- 

In general, for the unsaturated regime, it can be con- 
sidered that the rotation-age connection follows a Zipf-like 
power law ( Zipfj[T949, |1965| ), described as 

«(t) 

Vo 



1+'- 

T 



(10) 



where r is the timescale that controls the crossover between 
the initial plateau and the power-law domain, characterized 
by the exponent a, and vq is a normalization constant. Equa- 



tion ( 10 1 was deduced within Tsallis' nonextensive theory 
I Tsallis, Bemski fc Mendes|1999| [Lyra et al.\2003\ , exhibit- 
ing fractal properties and long-range memory. This theory 
demonstrates that the usual exponential for saturated fields 
can be recovered if a — > cx). 



j{t) = Vo exp [-Ai(i - to)] , 



(12) 



where Ai denotes the classical Lyapunov coefficient. 

According to Lyra et al- (20031, in contrast to the ex- 
ponential behavior presented in eq. (12 1, eq. (Ill must be 



replaced with a slower power law at critical points where 
long correlations develop. In this context, a similar expres- 
sion can be created to characterize the possible behavior of 
stellar rotation distribution in an unsaturated regime, de- 
noted by 



d_ 
dt 



Vsat 



-\a 



Vsat 



(A, ^0;q^ 1). 



Integrating the equation above, we have 



vit) 



[l + iq-l)\{t-tsat)]^-' 



(13) 



(14) 



where exp^a;) 



fl 



[1 + (1 - q)x 



1/(1-9) 



denoted by q- 



exponentiall and A, the nonextensive Lyapunov coefficient. 
Thus, in the nonextensive scenario, eq. ([7| is well de- 
scribed as a non-linear differential equation in form ( 13 1 with 
solution (14 1. When q — 1, eq. ( 12 I is recovered, indicating 



that a system governed by a saturated regime is in thermo- 
dynamic equilibrium. By contrast, when q differs from unity, 
a system controlled by unsaturation is out of equilibrium. 

In agreement with our model, for the unsaturated 
regime, exponent q can be described by a pair (a, A^) given 
by 



f 4- 



AaN 



(15) 



whereas for saturated regime g = 1 due to the usual expo- 
nential function. From above equations, we can estimate a 
relationship between the characteristic time ri and Tq given 
by 



T„ Al 



(16) 



in this context. As denote the braking strength. 

The g-parameter is related to the degree of nonextensiv- 
ity that emerges within thermodynamic formalism proposed 
in [Tsallis] ( |1988[ ). As revealed in eq. ( |15[ ), g is a function of 
magnetic field topology A'' and dynamo law a, which depend 



on stellar evolution. According to Kawaler (19881, small A'^ 



values result in a weak wind that acts on the MS phase of 
evolution, while winds with large A'^ values remove signifi- 
cant amounts of angular momentum early in the Pre-Main- 
Sequence (PMS). In this phase, for a given mass, maximum 
rotational velocity v decreases as A*' increases. This result is 
significant because, within a thermostatistical framework in 
which eq. ( |13[ ) naturally emerges, for a given value (a, N) we 
obtain the scale laws found in the literature, such as those 



proposed by Skumanich ( 1972 1 and Soderblom et al. ( 1991 1 



3.1 Saturated and unsaturated regimes 



Applying Tsallis statistics (Tsallis 20041, we can assume 
that eq. (|4| follows a simple linear differential equation in 
the form 



di \vo 
from which 



= -Ai 



(11) 
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3.2 /^-Saturated regime 



As reported by |Ivanova and Taam| ( [2003[ ), high rotation rates 
can differ significantly from the exponential decay, without 
which the angular momentum loss rate necessarily saturates 
the magnetic field (Chaboye r et al.\1995\ - As such, we used 
a more general relationship than that described by eq. (II]). 



See Borges 1 2004 1 for a review of g-calculus. 
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Figure 1. Semi-log plot of the rotational velocity distribuition for the sample of F and G-type stars. Solid lines are refered to usini 
data, while short dashed lines are the data multiplied by the statistical correction factor 4/n. 



Thus, this equation can be rewritten considering that the 
dJ/di rate depends on the SI, as follows: 



— = -K^.n'^uj^^t 



M 



-N/3 



M 
10-1' 



/3 l+4,aN/3-fl I R 

\Rq, 

1-2JV/3 

, {to^t< tsat), (17) 



for /3 = 1, we have the saturated domain, while for /3 > 1 
can be describe as a quastsaturated regime, defined here as 
/3-saturated. Deriving eq. ([3| and comparing with eq. ( |17[ ), 
we obtain that 



v{t) = Vo 
with 



it ~ to) 



T-/3 



1-/3 



(18) 



Tfj 



(/3 - 1)/(A) 



, /3-1 
yo \ 4aN/3 

Vsat 



Following the same procedure used to obtain v{t) in the 
previous cases, we have 



V{t) = wo [1 + (/? - l)Xp{t - to)] T^ , 



(20) 



where, /? is defined in interval 1 ^ /? < 1 + AaN/S. 

In this regime, a generalized expression to H can 
be obtained using the definition of so-called g-logarithm: 
ln,(x) = 2^^. Thus, 

tsat=to + {P~l)rfilnJ^^), (21) 



where for 13 = 1 the standard logarithm is recovered in ([9|. 
In this general case, we can extrated the rate between 
the time r, and r^, where we obtain 



TH 



q-l\ \ 



A/3 



(19) !^ 



where n and rp are related by 

1-/3 






P- 



vo 

Vsat 



(22) 



(23) 



if /3 = 1 , we recovered the rate between ri and Tq in case of 
the unsaturated-saturated model. 
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Figure 2. Semi-log plot of the stellar rotation-age relationship for the sample of F and G-type stars used in the present study. Upper 
and lower panels display mean projected rotational velocity, (vsmi), by age intervals of 0.5 Gyr for F and G-type stars, respectively, 
represented by the solid circles. In both panels, continuous lines represent the fit using Skumanich's law, and dash-dotted and dashed 
lines correspond to the fits using eqs. (|12[l and l|14|, respectively. Error bars for vsini and age are also indicated. 



4 THE STELLAR SAMPLE AND 
OBSERVATIONAL DATA 

Our study is based on the kinematically unbiased and 
magnitude-limited sample of nearby F and G dwarf stars 
from the Geneva-Copenhagen Survey (CGS) of the Solar 
neighbourhood, volume complete to ~ 40 pc, carried out 



by Nordstrom et al. ( 2004 1 and rediscussed by Holmberg, 



Nordstrom & Andersen (2007), which contain age, metallic 



ity, mass, projected rotational velocity v sini and kinematic 
properties for about 14000 F and G dwarfs. Readers are re- 
ferred to these authors for observational procedure and data 
reduction, but let us briefly discuss the quality of projected 
rotational velocity v sinj, age and mass measurements, given 
that these parameters are the most sensitive in the present 
study. 

For the vast majority of stars (12941 stars) the v sini 
measurements were computed from observations carried out 
with the CORAVEL spectrometers! Baranne et al. 1979 



of 1 km s , at least for rotat ions lower than about 30 km 
s"^ ( |De Medeiros fc Mayor|l999J . 

A total of 833 stars rotation measurements were deter- 
mined from observations carried out with the digital spec- 
trometer ( |Lathani||1985 l at the Harvard-Smithsonian Cen- 
ter for Astrophysics, CfA, in particular for stars rotating 
too rapidly for CORAVEL. For these stars, the v sini val- 
ues were computed on the basis of the best-fitting template 
spectrum ( [Nordstrom et aL|1994[ |1997[ ) . As shown by these 
authors, for slowly-rotating stars, typically for v sini values 
lower than 20 km s~^, there is a very good agreement be- 
tween data based on CfA and CORAVEL observations. As 
the vast majority of the present stellar working sample have 
rotation below 20 km s~^, we can consider that data on ro- 
tation used in the present analysis have the same quality, 
with a precision of about 1 km s"'^. 

We assumed that the mean value for a random distri- 



bution of inclinations sin i is given by 4/7r ( Chandrasekhar 



& Munch 19501. As mentioned by Reiners and Mohanty 



Mayor|1985l, applying (Benz fc Mayor_1980 1984l's calibra- (20121, this correction is only necessary in a statistically av- 



tions. Such a procedure gives v sini values with a precision 



eraged sense and not for comparison of single observed pop- 
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Figure 3. The Effective Temperature dependences of q (upper panel), A, (lower panel) used to fit the our sample segregated by different 
interval of spectral type. Upper figure indicates a linear dependende between the effective temperature and parameter l/(g — 1). For 
instance, Schumanich law (l/(q — 1) = 0.5) corresponds to stars with solar effective temperature. While that Pace & Pasquini law 
corresponds to hottest stars. 



ulations. Such a fact is valid for the present stellar sample, 
which is unbiased and magnitude-complete. In this sense 
statistical variations in the sin i distribution should not sig- 
nificantly change the presence or location of the very steep 
transition from a large population of undetected v sin i to a 



similarly large population of high v sin i observed ( Reiners 
and Mohanty|20T2|). Fig.Illshows the projected velocity dis 



tribution per bin uncorrected and corrected by factor 4/7r, 
in which no significant sin i effects can be observed. 

Good age measurements are, of course, crucial when 
investigating the rotation-age relationship. Most of the pre- 
vious studies on this subject have used ages estimated from 
chromospheric activity. Nevertheless, such a procedure is not 
valid for stars at the age of the Sun, because chromospheric 
emission essentially vanishes at this stage. The theoretical 
isochrone age is then the best choice, although important 
discrepancies subsist, in particular for old low-mass stars. 
The ages given in the GCS were then obtained on the basis 



of the stellar isochrone procedure Nordstrom et al. ( 2004 



using the Bayesian computational technique of _j0rgensen 
and Lindegren (20051. A revision of the data was carried 



sen| 



out by Holmberg, Nordstrom & Andersen ( 2007 1, taking into 



account new temperature and metallicity calibrations. Ac- 



cording Nordstrom et al. (20041 and Holmberg, Nordstrom 



& Andersen (20071 for 81% of the presumably single stars 



in the original sample of the CGS, age estimations based, 
on their procedure, have errors below 50%. 



Nordstrom et al. { 2004 1 have also estimated stellar 



masses using theoretical isochrone analysis, with individ- 
ual errors averaging about 0.05 Mq, and metallicities from 
Stromgren ubvyB photometry. As shown by these authors, 
the distribution of their estimated metallicities obeys a 
Gaussian distribution with a mean of -0.14 and disper- 
sion 0.19 dex. The distribution of the metallicities from 



IHolmberg, Nordstrom & Andersen ( 2007 1 for the 5835 sin 



gle stars used in the present analysis shows a mean value 
{[Fe/H]} = —0.18 and dispersion of 0.21 dex. Clearly, one 
observes that our working sample is mostly composed of 
solar-type stars, in spite of the presence of a low percent- 
age of stars with slightly sub-solar metallicities. Recently, 
Casagrande et al. (2011 1 re-analyzed most of the parameters 



for the CGS via the infrared flux method (IRFM). Accord- 



© 2002 RAS, MNRAS 000, [T]-?? 



A nonextensive approach to the stellar rotational evolution I 7 



ing to these authors, this method, in comparison to previ- 
ous CGS cahbrations, reveals differences in measurement of 
effective temperature leading to a much better agreement 
between stars and isochrones in the HR diagram. 

In order to test our nonextensive approach for stellar 
rotational magnetic breaking, we used v sini and ages from 
the referred Geneva-Copenhagen survey for F- and G-type 



single stars (Nordstrom et al. 2004 Holmberg, Nordstrom 
fc Andersen|2007l, with ages limited to 10 Gyr and masses 



of O.9OM0 sj M ^ 2.OM0. This upper limit for ages was 
defined by the star's lifetime on the MS and the age of the 
Galactic disk ( Epstein fc Pinsonneault|2012 |, while the up- 
per limit for masses (or spectral type) is based on the differ- 
ent internal structure of the stars in our sample. These age 
and mass limits were chosen to avoid the presence of stars 
exhibiting the most significant uncertainties for these pa- 
rameters, namely the oldest and lowest mass stars, as well 
as the highest mass stars. Based on these criteria, we ob- 
tained a final working sample of 5780 single stars consisting 
of 3790 F- (2050 with spectral type between FO and F5 and 
1750 between F6 and F9) and 1990 G- (1919 with spectral 
type between GO and G5 and 73 between G6 and G9) type 
stars. 



5 TESTING THE MODEL WITH 
UNSATURATED STARS 

First of all, stars in the present working sample can be con- 
sidered to be in the unsaturated regime , indicating that the 
computation of the values of the parameter /3 is not nec- 
essary. Indeed, all the stars are aged more than 1 Gyr, as 
described in Section 4. 



Different studies in the literature {e.g: ( Krishnamurthi 



et a/.|1997||Bouvier, Forestini fc Allain|1997| [Epstein fc Pin 



sonneault|2012 1) use the rotation distribution of young open 
clusters to estimate initial conditions. Here, we follow the 
same approach, by considering the rotation distribution of 
T Tauri stars. In agreement with initial conditions, we com- 
pute the best-fit [g,Aq] doublet. Fig. [2] and [3] shows the re- 
sults of the present analysis, where the mean of v sin i by age 
intervals of 0.5 Gyr are shown for the aforementioned sample 
of F and G-type stars. In Fig. [2J the upper and lower pan- 
els also display the best fits applying the Skumanich law, as 
well as eq. ( |14[ ). The best-fit parameters for the rotation-age 
relationship using this equation are depicted in Fig. |3] 

Fig. 12] clearly shows that the saturated model does not 
fit the data for either type of star aged > 1 Gyr. This sug- 
gests that the stars in our sample are still under the effect 
of magnetic torque. Additionally, there is a positive gradi- 
ent for index q from the hottest to the coolest stars, namely 
from F- to G-type stars, suggesting that at a given age the 
coolest stars have a stronger rotational decline, reaching the 
slow rotation regime faster than the hottest stars. More- 
over, there is also a positive gradient for braking strength 
Ag from the hottest to the coolest stars. Indeed, from Table 
[1] demonstrates that the present model clearly depend on 
stellar mass. 

According to q values presented in Fig. [S] magnetic 
field geometry is essentially (quasi) dipolar, irrespective of 
the value of a. For this type of topology, our results indi- 
cate that for F early-type stars the Pace and Pasquini law 



{q = 1.4) is the most appropriate, while for late- type stars 
the Soderblom et al.'s {q — 2.3) law can be considered a good 
approach. In case of G stars the Skumanich {q — 3) law is 
the most appropriate for early-type and the Soderblom et 
al. law is for late-type. 



6 CONCLUSIONS 

In summary, our study presents a new statistical approach 
for stellar rotational braking. The present analysis shows 
that the rotation-age relationship can be well reproduced 
using a nonextensive approach from statistical mechanics, 
namely the Tsallis nonextensive models. According to dif- 
ferent studies {e.g 



Rutten fc Pylyser (19881; Matt et al. 



(20121), a theoretical model for magnetic braking should 



consider several parameters, {e.g. magnetic field geometry, 
mass loss rate and coronal temperature), which to data are 
not represented in the different laws describing the rotation- 
age relations. The nonextensive analysis proposed here con- 
siders some of these parameters, and offers the possibility of 
studying the stellar rotational braking behaviur of different 
classes of stars, based on the same approach. 

Our theoretical-observational results reveal that each 
spectral type range is characterized by a specific entropy 
production, this suggests the occurrence of different rota- 
tional velocity probability distribution with high energy tails 



as proposed by de Freitas & Be Medeiros (20121. Moreover, 



these results suggest that increasing index q behavoir be- 
tween the saturated and unsaturated regimes as a function 
of spectral type may be due to a phase transition magnetic 
process that occurs in the stellar dynamo dynamics. Accord- 



ing to Bouvier, Forestini & AUain (19971, low mass stars 



retain the memory of their initial angular momentum up to 
an age of about 1 Gyr, as opposed to solar type stars, where 
the memory is retained up to few 10** years. This memory 
can be quantified using the entropic index q. The fact that 
this index increases from F to G stars, indicates that long- 
range memory grows at the same scale. Higher mass stars do 
not maintain this type of memory for long time, because the 
interaction of the young star with its circumstellar disk is 
much shorter than their very low-mass counterparts, which 
retain constant angular velocity. 

Finally, a comparison with observational data, namely 
projected rotational velocity wsini, shows that our nonex- 
tensive approach closely fits the decrease in stellar rotation 
with age for F- and G-type main-sequence stars. Indeed, 
the behavior of the index q, which is related to the degree of 
nonextensivity, and in the present aproach can be described 
by a and A*', associated to the dynamo process and mag- 
netic field geometry, respectively offers relevant information 
concerning the level of magnetic braking. At a given age the 
coolest stars have a stronger rotational decline, reaching the 
slow rotation regime faster than the hottest stars. 
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